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We compute the partition function of the su(m) Polychronakos-Frahm spin chain of BCn type by 
means of the freezing trick. We use this partition function to study several statistical properties of 
the spectrum, which turn out to be analogous to those of other spin chains of Haldane-Shastry type. 
In particular, we find that when the number of particles is sufficiently large the level density follows 
a Gaussian distribution with great accuracy. We also show that the distribution of (normalized) 
spacings between consecutive levels is of neither Poisson nor Wigner type, but is qualitatively similar 
to that of the original Haldane-Shastry spin chain. This suggests that spin chains of Haldane- 
Shastry type are exceptional integrable models, since they do not satisfy a well-known conjecture of 
Berry and Tabor according to which the spacings distribution of a generic integrable system should 
be Poissonian. We derive a simple analytic expression for the cumulative spacings distribution of 
the BGjv-type Polychronakos-Frahm chain using only a few essential properties of its spectrum, like 
the Gaussian character of the level density and the fact the energy levels are equally spaced. This 
expression is in excellent agreement with the numerical data and, moreover, there is strong evidence 
that it can also be applied to the Haldane-Shastry and the Polychronakos-Frahm spin chains. 

PACS numbers: 75.10.Pq, 05.30.-d, 05.45.Mt 



I. INTRODUCTION 

Solvable spin chains often provide a natural setting for 
testing or modeling interesting physical phenomena and 
mathematical results in such disparate fields as fractional 
statistics, random matrix theory or orthogonal polynomi- 
als. Among these chains, those of Haldane-Shastry (HS) 
type occupy a distinguished position due to their remark- 
able integrability and solvability properties. The original 
chain of this type was independently introduced twenty 
years ago by Haldane [l[ and Shastry 0, in an attempt 
to construct a model whose ground state coincided with 
Gutzwiller's variational wave function for the Hubbard 
model in the limit of large on-site interaction III In 
the original HS chain, the spins are equally-spaced on a 
circle and present pairwise interactions inversely propor- 
tional to their chord distance. 

An essential feature of the spin chains of HS type 
is their close connection with the spin versions of the 
Calogero [|| and Sutherland models, and their gen- 
eralizations due to Olshanetsky and Perelomov This 
observation — already pointed out by Shastry in his orig- 
inal paper — was elegantly formulated by Polychronakos 
in Ref. In the latter reference, the author showed 

that the original HS chain can be obtained from the 
spin Sutherland model [ill El, EH m the strong cou- 
pling limit, in which the dynamical and spin degrees of 
freedom decouple, so that the particles "freeze" at the 
equilibrium positions of the scalar part of the potential. 
In this regime, the integrals of motion of the spin Suthcr- 
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land model directly yield first integrals of the HS chain, 
thereby explaining its complete integrability. This pro- 
cedure was also applied in [l(| to construct a new inte- 
grable spin chain related to the original Calogero model. 
The spectrum of this chain was numerically studied by 
Frahm (l4| . who found that the levels are grouped in 
highly degenerate multiplets. In a subsequent publica- 
tion, Polychronakos computed the partition function of 
this chain (usually referred to in the literature as the 
Polychronakos-Frahm chain) by the "freezing trick" ar- 
gument described above [1 51 ] - Interestingly, the partition 
function of the original HS chain was computed only very 
recently [lj|. 

Both the HS and the PF (Polychronakos-Frahm) 
chains are obtained from the Sutherland and Calogero 
models associated with the An root system in Olshanet- 
sky and Perelomov's approach. The BCn versions of 
both chains have also been studied in the literature. More 
precisely, the integrability of the PF chain of BCm type 
was established by Yamamoto and Tsuchiya [171 ] using 
again the freezing trick. On the other hand, the parti- 
tion function of the HS chain of BCn type was computed 
in closed form in Ref. [l8| . The explicit knowledge of the 
partition function made it possible to study certain sta- 
tistical properties of the spectrum of this chain. In par- 
ticular, it was observed that for a large number of spins 
the level density is Gaussian. As a matter of fact, this 
property also holds for the original HS chain, as shown 
in Ref. [16| . The analysis of the distribution of the spac- 
ing between consecutive levels of the original HS chain 
was also undertaken in the latter reference. Rather un- 
expectedly, it was found that this distribution is not of 
Poisson type, as should be the case for a "generic" inte- 
grable model according to a long-standing conjecture of 
Berry and Tabor This behavior has also been re- 

cently reported for a supersymmetric version of the HS 
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chain |20|. 

The aim of this paper is twofold. In the first place, 
we shall compute in closed form the partition function 
of the PF chain of BCn type by means of the freezing 
trick. Using the partition function, we shall perform a 
numerical study of the density of levels and the distri- 
bution of the spacing between consecutive energies. We 
shall see that the level density is again Gaussian, and 
that the spacings distribution is analogous to that of the 
original HS chain. In particular, our results show that 
the distribution of spacings is neither Poissonian nor of 
Wigner type (characteristic of chaotic systems). We shall 
next derive a simple analytic expression for the cumula- 
tive spacings distribution, which reproduces the numer- 
ical data with much greater accuracy than the empiric 
formula proposed in Ref. [u|. In fact, we have strong 
numerical evidence that the new expression can also be 
applied to the HS and PF chains of An type. In view 
of the Berry-Tabor conjecture, our results suggest that 
spin chains of HS type are exceptional among the class 
of integrable models. 



BCn Calogero model 



H sc = 



•5X+a(a-l)£ 



(4) 2 



(3) 



with b = (3a and a > 1/2. Moreover, it can be shown [22| 
that f j = \/^yi, where ?/j is the i-th zero of the general- 
ized Laguerre polynomial L N ~ X . From this fact, one can 
infer [23|] that for N ^> (3 the density of sites (normalized 
to unity) Pn{x) is given by the circular law 



Pn(x) 



1 



2ttN 



V8N- 



(4) 



Note that in this limit the sites' density is independent 
of (3 and is qualitatively similar to that of the PF chain 
of An type [14| . Integrating the previous equation, we 
obtain the implicit asymptotic relation 



Auk = t k \l%N-£, 2 k + 8iVarcsin (-4=) 



II. THE PARTITION FUNCTION OF THE PF CHAIN 
OF BC N TYPE 



valid also for N ^> (3. 

The spin chain (fT]) can be expressed in terms of the 
spin Calogero model of BCn type 



The Hamiltonian of the (antiferromagnetic) su(m) PF 
chain of BCn type is defined by 



K 



E 



1 + 5, 



1 + 5, 



eSi 



where the sums run from 1 to N 
unless otherwise stated), (3 > 0, e 



(1) 



(as always hereafter, 
= ±1, Sij is the op- 



erator which permutes the i-th and j-th spins, Si is the 
operator reversing the i-th spin, and Si- 



Si Sj Sij . 



Note 



that the spin operators Sij and Si can be expressed in 
terms of the fundamental su(m) spin generators Ju at 
the site k (with the normalization ti(J% J^) = ^5 ai ) as 



^ = -EC+«E 

i i^j 



s,. 



K) 2 



*E 



and its scalar reduction (O as 



a - 



Si 



K) 2 



b — eS., 



are 



r 2 (5) 



(6) 



where 7i e is obtained from Tt e replacing the chain sites £ 
by the particles' coordinates x. Since 



H' 



dl +a 2 U + Q{a). 



Sij — 

m 



a=l 



S t = 



The chain sites £j are the coordinates of the unique min- 
imum in C = {x | < x\ < ■ ■ ■ < xn} of the potential 



£/(x) = E 



1 



1 



K) 2 



r 

T 



(2) 



where " 



± Xj and r 2 = J2i x f- The existence of 
this minimum follows from the fact that U tends to +oo 
on the boundary of C and as r — > oo, and its uniqueness 
was established in Ref. [2l[ by expressing the potential 
U in terms of the logarithm of the ground state of the 



when the coupling constant a tends to infinity the parti- 
cles in the spin dynamical model concentrate at the 
coordinates of the minimum of the potential U, that is 
at the sites £j of the chain (fTJ). Thus, in the limit a — > oo 
the spin and dynamical degrees of freedom of the Hamil- 
tonian (O decouple, so that by Eq. © its eigenvalues are 
approximately given by 



Ef 



a£* 



a > 1 



(7) 



where Ef c and £j are two arbitrary eigenvalues of Hf c 
and H e , respectively. The asymptotic relation ([7]) imme- 
diately yields the following exact formula for the partition 
function Z e of the chain ([TJ) : 



Z e (T) = lim 



Z e (aT) 
Z sc (aT) 



(8) 
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where Z e and Z sc are the partition functions of H e and 
H sc , respectively. 

We shall next evaluate the partition function Z e of 
the chain ([1]) by computing the partition functions Z e 
and Z sc in Eq. ([5]). In order to determine the spectra of 
the corresponding Hamiltonians H e and H sc , following 
Ref. [16j | we introduce the auxiliary operator 



: (a-K ij ) + 



(9) 



where Kij and A^ are coordinate permutation and sign 
reversing operators, defined by 

{-Kij f){ x li ■ ■ ■ i x ii ■ ■ ■ ,Xj, . . . , Xtf) 

f \X\ ) • • ' 5 ? * * • 5 ) * * • 3 *^iV ) 3 

(K l f)(x 1 , . . . ,Xi, . . .,x N ) = f(xx, . . .,-Xi, . . . ,x N ) , 



and = KiKjKij. We then have the obvious relations 



H e = H'\ 



H'\ 



(10a) 
(10b) 



On the other hand, the spectrum of H 1 can be easily 
computed by noting that this operator can be written in 
terms of the rational Dunkl operators of BCn type [24| 



Ji =d Xi + 



-4(1-^) + 4(1-^-; 



+ - (1 - Ki) , 

Xi 

as follows |25j : 

h' = [ - ( J r) 2 + a J2 Xid ** + E < 

i i 

where 

Mx)= e -t r2 ni^rn 



AT. 
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is the ground state of the Hamiltonian (|3|) and 
£ =iVa(Wa(iV-l) + -). 



(11) 



(12) 



(13) 
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Since the Dunkl operators (|11[) map any monomial Y\ i a;™* 
into a polynomial of total degree n% + ■ ■ ■ + npf — 1, by 
Eq. (fT2"j) the operator H' is represented by an upper tri- 
angular matrix in the (non-orthonormal) basis with ele- 
ments 



n 



n=(m,...,n N ) e (EU{0}f, (15) 



ordered according to the total degree |n| = n\ 
of the monomial part. More precisely, 



H'4> n = E' n (/} n 



|m|<|n| 



where 



E'=a\n\+E Q 



-n N 



(16) 



(17) 



and the coefficients c mn are real constants. 

We shall now construct a basis of the Hilbert space of 
the Hamiltonian H f in which this operator is also rep- 
resented by an upper triangular matrix. To this end, 
let us denote by A e the projector on states antisymmet- 
ric under simultaneous permutations of spatial and spin 
coordinates, and with parity e under sign reversals of co- 
ordinates and spins. If 



— m — 1 

— 2 1 



|s) = . . . , s N ), s t = -M, —M + 1, . . . , M 
denotes a state of the su(m) spin basis, the functions 



^ s (x)=A e (0 n (x)|s)), 



(18) 



form a basis of the Hilbert space of the Hamiltonian H e 
provided that: 

(i) Ui ^ • • • ^ npf. 

(ii) Si > Sj whenever m = rij and i < j. 

(iii) s l ^ for all i, and s t > if (-1)" ! = -e. 

The first two conditions are a consequence of the anti- 
symmetry of the states (|18p under particle permutations, 
while the last condition is due to the fact that these 
states must have parity e under sign reversals. Since 
KijA* = -SijA e and KiA e = eS l A e , it follows that 
H e A e = H'A e . Using this identity and the fact that 
H' obviously commutes with A £ , from Eq. (|16[) we easily 
obtain 

ff e VVs = #Vn,B = A £ ((#> n )|s» 

= E' n 1p ntS + ^ c mnV'm, S ■ 
|m|<|n| 

Thus H e is represented by an upper triangular matrix in 
the basis (fT5)) . ordered according to the degree |n|. The 
diagonal elements of this matrix are given by 



E n 



(19) 



where n and s satisfy conditions (i)-(iii) above. Note 
that, although the numerical value of E^ s is independent 
of s, the degeneracy of each level clearly depends on the 
spin through the latter conditions. 

Turning next to the scalar Hamiltonian H sc , in view 
of Eq. (|10b[) we now need to consider scalar functions of 
the form 



V> n (x) = A s ^ n (x) , 



(20) 



where A s is the symmetrizer with respect to both per- 
mutations and sign reversals. These functions form a 
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(non-orthonormal) basis of the Hilbert space of H sc pro- 
vided that rii — 2ki are even integers and k± ^ • • • ^ kpf. 
Just as before, the matrix of the scalar Hamiltonian H sc 
in the basis (|20[) ordered by the degree is upper triangu- 
lar, with diagonal elements E^ also given by the RHS 

of dm). 

Let us next compute the partition functions Z sc and 
Z e of the models ([3]) and ([5]). To begin with, from now 
on we shall drop the common ground state energy E 
in both models, since by Eq. ([5]) it does not contribute 
to the partition function Z e . With this convention, the 
partition function of the scalar Hamiltonian _ff sc is given 
by 



Z sc {aT) 



E n 



2|k| 



where q — e 1 /( k & T ). The previous sum can be evaluated 
by expressing it in terms of the differences pi — ki — fci+i, 



N 



i = 1, . . . , N — 1, with pn — fcjv- Since kj — Pi, we 



easily obtain 
Z sc (aT) = 



N N 

2 E Ep. 

j = l i=j 



E « 

pi ,...,pw>0 

n e (<?t - - ^ 



E 9 i=l 

pi,...,pw^0 
2i\-l 



(21) 



i pt^O 



In order to compute the partition function of the spin 
Hamiltonian H e , we shall first assume that m is even, so 
that condition (iii) simplifies to 

(hi') Si > for all i. 

As neither the value of E^ s nor conditions (i), (ii) and 
(iii') depend on e, in this case the partition functions Z e 
and Z e cannot depend on e. Hence from now on we shall 
drop the superscript e when m is even, writing simply 
Z and Z. By Eq. (p~9|) . after dropping E the partition 
function of the Hamiltonian ([5]) can be written as 



Z(aT)= d »9 |n1 ' ( 22 ) 

ni > ■ ■ -^n^v 5^0 



where the spin degeneracy factor d n is the number of spin 
states |s) satisfying conditions (ii) and (iii'). Writing 



n = ( ki, . . . , ki, . . . ,k r , ■ ■ ■ , k r ), k\ > ■ ■ ■ > k r ^ 0, 
by conditions (ii) and (iii') we have 



*. -n ("£)-«■" 

i— 1 x 7 



v = (ui, ...,u r ). (23) 



Note that ^1=1 = so that the multiindex v can be 
regarded as an element of the set Vn of partitions of N 



(taking order into account). With the previous notation, 
Eq. (HH) becomes 

Z(aT)= <*(") E 

vEVn ki>— >k r ^0 

= E d M E ^ 31 

vev N pi,...,p,._i>o 



e^) e n 

v€Pn pi,...,p r _i>0 j'=l 

Pr^O 



Pj E 



^E^)n T 

t^ePiv 3=1 



(24) 



where 



Nj = j> . 

i=i 

From Eqs. ©, (J2TJ) and $Q we finally obtain the follow- 
ing explicit expression for the partition function of the 
su(m) PF chain of BCpj type in the case of even m: 



<-V) at. 

Z{T) = q~ N Y[(l - q 2i ) d ( u ) II i~T^7 ' ^ 

i vev N j=i 

where t(y) = r is the number of components of the mul- 
tiindex v . For instance, for spin 1/2 we have Vi = 1 for 
all i, and therefore r — N, d(i/) = 1 and iVj — j, so that 
the previous formula simplifies to 

Z(T) = g^(JV-i) J](l + q *) , m = 2. (26) 



Thus, for spin 1/2 the spectrum of the chain {!]) is given 
by 

£j = ~ N(N-l)+j , j = 0, 1, .. \ N(N+1) , (27) 

and the degeneracy of the energy £j is the number Qn{j) 
of partitions of the integer j into distinct parts no larger 
than N (with Q N (0) = 1). For j sC TV this number 
coincides with the number Q(j) of partitions of j into 
distinct parts, which has been extensively studied in the 
mathematical literature 26]. It is also interesting to ob- 
serve that the partition function (|2"6"]) is closely related to 
Ramanujan's fifth order mock theta function (2?| 



N=0 



where Zm (q) denotes the RHS of Eq. 

Equation (f2"6"|) shows that for spin 1/2 the chain ((T|) 
is equivalent to a system of N species of noninteracting 
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fermions (with vacuum energy £q = N(N — l)/2), whose 
effective Hamiltonian is given by 



N 



H, 



cff — Ml 



i=l 



Here 6, (resp. b\) is the annihiiation (resp. creation) op- 
erator of the i-th species of fermion, and Ei = i its 
energy. A similar result was obtained in Ref. [28[ for 
the supersymmetric su(l|l) (ferromagnetic) HS chain, al- 
though in the latter case the energy of the i-th. fermion 
is Ei = i(N — i) (the dispersion relation of the original 
Haldane-Shastry chain) . 

Let us consider now the case of odd m. In this case, it is 
convenient to slightly modify condition (i) above by first 
grouping the components of n with the same parity and 
then ordering separately the even and odd components. 
In other words, we shall write n = (n c , n Q ), where 



n c = (2fci,...,2fci,.. 

V B +1 



, 2k s , . . . , 2k s J , 



; 2k s+1 + 1, . . . , 2k s+1 + 1, . . . , 2A; r + 1, . . . , 2k r + 1), 



and 



k\ > ■ ■ ■ > k s ^ 0, k s+1 > ■ ■ ■ > k r ^ . 
By conditions (ii) and (hi), the spin degeneracy factor is 



d n=n 



1=1 



2 



n 

:=s+l 



m — € 
2 \ — 



(28) 



Calling 



*i= E 



3 = s + 1, . . . ,r , 



i=s+l 

and proceeding as before, we obtain 



Z e (aT)= J2 E d » E ^ 



k 1 >->k s ^0 
fc s+ i>->fe r >0 



E E'»« 



N r 



E 9 

fel>--->fes>0 



E 2i/ i fc i 



E 



fe S+ l> -->feT-^0 



l{y) 

E E'-M' 

we^K s=o 



,-(N+N.) 



E 2 ^i fc i 



JIN; 



4=1 



n.r 



-,2Ni 



(29) 



J = S + 1 



Substituting the previous expression and (j2Tj) into (5]), 
we immediately deduce the following explicit formula for 
the partition functions of the su(m) PF chain of BCm 
type for odd to: 



2 e (T) = na 



I 21 ) E Y. d s(»)i- {N+Ns) 



4=1 



n.r 



9 



2Wi 



(30) 



Although we have chosen, for definiteness, to study 
the antiferromagnetic chain ([T]) , a similar analysis can be 
performed for its ferromagnetic counterpart 



T~tv — 



E 

»#4 



1 



1 - Si 



i 



C.S; 



Since now 



H f = H '\ K 



Su.Ki-mSi 



(31) 



(32) 



we must replace the operator A £ in Eq. (|18|) by the pro- 
jector on states symmetric under simultaneous permuta- 
tions of the particles' spatial and spin coordinates, and 
with parity e under sign reversal of coordinates and spin. 
Hence condition (ii) above on the basis states ip n ,s should 
now read 

(ii') Si Sj whenever m = Uj and i < j. 

As a result, the degeneracy factors d(v) and d e (v) in 
Eqs. (p?3"|) and (|28p should be replaced by their "bosonic" 



s 

*w = n 



Vi - 1 



^ + Vi - 1 



n 

i=s+l 



+ !/,--! 



Therefore the partition function of the ferromagnetic 
su(m) PF chain of BC N type (J3TJ) is still given by Eq. $25j) 
(for even m) or ([30]) (for odd to), but with d(v) and d e (f ) 
replaced respectively by d-p{u) and dp^). 

On the other hand, the chains |T]) and (|31 [) are obvi- 
ously related by 

^+^- £ = 2(j2 [te-e,r 2 +te+^r 2 ]+/3E^ 

The RHS of this equation clearly coincides with the 
largest eigenvalue fmax of the antiferromagnetic chains 
7i £ , whose corresponding eigenvectors are the spin states 
symmetric under permutations and with parity e under 
spin reversal. This eigenvalue is most easily computed 
for the spin 1/2 chains, since in this case the spectrum is 
explicitly given in Eq. (f2"T|) . We thus obtain 

£max = I N(N + l N(N + 1)=N\ (33) 
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so that 

Tip = N 2 — . 

Hence the partition functions Z e and Zp of H e and Til- 
satisfy the remarkable identity 

Z^( q )=q N2 Z^( q - 1 ). 

This is a manifestation of the boson-fermion duality dis- 
cussed in detail in Ref. [29( for the su(m|n) supersym- 
metric HS spin chain, since the ferromagnetic (resp. an- 
tiferromagnetic) chain can be regarded as purely bosonic 
(resp. fermionic). For instance, using the latter identity 
and Eq. (|2"6")l we easily obtain the following expression 
for the partition function of the ferromagnetic spin 1/2 
chains: 

ZF(T)=l[(l + q l ), m = 2. (34) 

i 

(Note that, as in the antiferromagnetic case, Z F is actu- 
ally independent of e for even m.) This is, again, the par- 
tition function of a system of N species of free fermions 
of energy Ei = i, but now the vacuum energy vanishes. 

Equation (|25| for the partition function of the antifer- 
romagnetic chains with even m can be easily simplified 
to 

Z(T) = + <f ) £ d(y) q '=* U (l~q N *), 
i vev N 3 = 1 

where the positive integers N'- are defined by 

{N{, N' N _ eH } = {1, . . . , N-l}-{Ni, N l(v) ^}. 

The sum in the RHS is easily recognized as the parti- 
tion function Z l > A * > (T;m/2) of the su(m/2) (antiferro- 
magnetic) PF chain of An type [2(J. We thus obtain 
the remarkable factorization 

Z(T;m) = Z F (T;2) • Z (A) (T; m/2) , m 6 2N , (35) 

where the second argument in Z and Z-p denotes the 
number of internal degrees of freedom. Replacing d(v) 
by dp{y) in Eq. (|25p we obtain a similar factorization for 
the partition function of the ferromagnetic chains: 

Z F (T;m) = Z F (T;2) ■ Z F A) (T; m/2) , m e 2N . (36) 

Thus, for even m the PF chains of BCm type |T]) and (|3T|) 
can be described by an effective model of two simpler 
nonintcracting chains. This remarkable property, which 
to the best of our knowledge is unique among the class of 
chains of Haldane-Shastry type, certainly deserves fur- 
ther investigation. 



III. SPACINGS DISTRIBUTION AND THE 
BERRY-TABOR CONJECTURE 

For fixed values of the number of particles N and the 
internal degrees of freedom m, it is straightforward to 
obtain the spectrum of the chain ([T]) by expanding in 
powers of q the expressions (|25|) or (|30|) for its partition 
function. In this way, we have been able to compute the 
spectrum of the latter chain for relatively large values of 
N (for instance, up to N = 22 for m — 3). Our calcu- 
lations conclusively show that the spectrum consists of 
a set of consecutive integers. For even m, this observa- 
tion follows immediately from the expression (|34[) - (|3"5]) 
and the fact the the energies of the PF chain of Am type 
are also consecutive integers. For odd m we have been 
unable to deduce this property from Eq. (|3"0|) for the par- 
tition function, although we have verified it numerically 
for many different values of N and m. 

Our computations also evidence that for TV > 10 the 
level density (normalized to unity) can be approximated 
with great accuracy by a normal distribution 

1 (£-(i) 2 

g(S) = - r e-^-, (37) 

where fj, and a are respectively the mean and the variance 
of the energy. For instance, in Fig. [1] we compare the 
cumulative level density 

F(£) = m- N dt, 



where Ei is the i-th energy and di its degeneracy, with 
the cumulative Gaussian density 




140 160 180 200 220 



FIG. 1: Cumulative distribution functions F(£) (at its dis- 
continuity points) and G(£) (continuous red line) for TV = 15 
and spin 1/2. 
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the approximately Gaussian character of the level density 
has already been verified for other chains of HS type, 
like the original Haldane-Shastry chain [l6j], its su(m|n) 
supersymmetric version [20| , and the HS spin chain of 
BC N type 0. 

The mean energy fi and its standard deviation a 2 , 
which by the previous discussion characterize the approx- 
imate level density of the chain ([T]) for large N, can be 
computed in closed form. Indeed, in Appendix [X] we 
show that 



l(l + IW-^-(l + ep), 
2 \ ml Im 



N 



_(4A 2 + 6A-l)(l- — 



1 



N(l-p) 
Am 2 



(39) 
, (40) 



where p S {0, 1} is the parity of m. Thus, when N tends 
to infinity \i and a 2 respectively diverge as N 2 and iV 3 , 
as for the original Polychronakos-Frahm chain [3|| . By 
contrast, it is known that \i ~ TV 3 and a 2 ~ A*" 5 for the 
trigonometric HS chains of both An [n| and BCn [IH 
types. It is also interesting to observe that the standard 
deviation of the energy is independent of e even for odd 
m, when the spectrum does depend on e according to the 
previous section's results on the partition function. 

We have next studied the probability density p(s) of 
the spacing s between consecutive (unfolded) levels of 
the chain (TTJ). For many important integrable systems it 
is known that p(s) is Poissonian [33, f3lj ] . in agreement 
with a well-known conjecture of Berry and Tabor [l9j |. 
On the other hand, it has been recently shown that for 
the HS chain of An type [l6| (and its supersymmetric 
extension 20]) the cumulative density P(s) = p(x)dx 
is well approximated by an empiric law of the form 



P(s) 



1-7 1- 



(41) 



where s max is the largest normalized spacing and a, (3 are 
adjustable parameters in the interval (0, 1). The param- 
eter 7 is fixed by requiring that the average spacing be 
equal to 1, with the result 
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(— + ^ + (42) 



where B is Euler's Beta function. Thus, the cumulative 
density of spacings for the HS chains of A n type follows 
neither Poisson's nor Wigner's law 

P( s ) = i-e-™ 2 / 4 , 

characteristic of a chaotic system. Our aim is to ascertain 
whether the cumulative density of spacings for the PF 
chain of BCn type JT]) resembles that of the Ajy-type HS 
chain, or is rather Poissonian as expected for a generic 
integrable model. 

In order to compare the spacings distributions of spec- 
tra with different level densities, it is necessary to trans- 
form the "raw" spectrum by applying what is known as 



the unfolding mapping [32| . This mapping is defined by 
decomposing the cumulative level density F{£) as the 
sum of a fluctuating part Ff\{£) and a continuous part 
t)(£), which is then used to transform each energy Si, 
i = l,...,n, into an unfolded energy rji — T)(£i). In 
this way one obtains a uniformly distributed spectrum 
{Vi}2=ii regardless of the initial level density. One finally 
considers the normalized spacings Sj = (rji+i — rji) /A, 
where A = (r) n — 771) / (n — 1) is the mean spacing of the 
unfolded energies, so that {sj}^ has unit mean. 

By the above discussion, in our case we can take the un- 
folding mapping rj{£) as the cumulative Gaussian distri- 
bution (j3"8|) . with parameters ii and a respectively given 
by ([39]) and (|40|) . Just as for the level density, in order 
to compare the discrete distribution function p(s) with 
a continuous distribution it is more convenient to work 
with the cumulative spacings distribution P(s). Our 
computations show that for a wide range of values of N, 
m and e = ±1, the distribution P(s) is well approximated 
by the empiric law (|4"Tj) with suitable values of a and f3. 
For instance, for N — 20 and m = 2 the largest spacing is 
•Smax = 3.13, and the least-squares fit parameters a and /3 
are respectively 0.29 and 0.24, with a mean square error 
of6.0xl0~ 4 (see Fig. g). Thus the PF spin chain of BCn 




FIG. 2: Cumulative spacings distribution P(s) and its ap- 
proximation P(s) (continuous red line) for iV = 20 and m = 2. 
For convenience, we have also represented Poisson's (green, 
long dashes) and Wigner's (green, short dashes) cumulative 
distributions. 

type behaves in this respect as the HS chain of An type, 
and unlike most known integrable systems. In fact, we 
have also studied the spacings' distribution of the orig- 
inal (Ajy-type) PF chain, obtaining completely similar 
results [35]. These results (and also those of Ref. [2pj ) 
suggest that a spacings distribution qualitatively similar 
to the empiric law (|4"Tj) is characteristic of all spin chains 
of HS type. 

Our next objective is to explain this characteristic be- 
havior of the cumulative spacings distribution P(s) of the 
chain |T]) using only a few essential properties of its spec- 
trum. We shall find a simple analytic expression without 
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any adjustable parameters approximating P(s) even bet- 
ter than the empiric law (flTj) . Moreover, we have strong 
numerical evidence that the new expression also provides 
a very accurate approximation to the cumulative spacings 
distribution of the original HS and PF chains. 

Consider, to begin with, any spectrum £\ = £ m in < 
■■■<£„ = £ max obeying the following conditions: 

(i) The energies are equispaced, i.e., — = 5 for 
i = 1,.. . ,n — 1. 

(ii) The level density (normalized to unity) is approxi- 
mately Gaussian, cf. Eq. ([57)1 . 

(iii) £ n - (i , il - £ i > a. 

(iv) £ i and £ n are approximately symmetric with respect 
to /i, namely \£i + £ n — 2fi\ <C £ n — £±. 

As discussed above, the spectrum of the chain {T]) satisfies 
the first condition with 5 = 1, while condition (ii) holds 
for sufficiently large N. As to the third condition, from 
Eqs. (03]), (0S|), (HI, jEEJ and jB2| it follows that both 
(£ n ~ fi)/cr and (fj, — £i)/cr grow as N 1 / 2 when N —> oo. 
The last condition is also satisfied for large N, since by 
the equations just quoted \£i + £ n — 2/i| = O(N) while 
£ n -£ x = 0(N 2 ). 

From conditions (i) and (ii) it follows that 

- m = G(£ l+1 ) - G(£i) ~ G'(£i)S = g(£i)S. 

On the other hand, by condition (iii) we have 

J7« = G(£„) ~ 1 , ryi = G(5i)~0, 

so that A = l/(n - 1). Thus 

Vi+i - Vi 



Si 



where 



Wg{£i) 



e I*' 2 



2na 



W=(n-l)5 = £ n - £x 



(43) 



(44) 



on account of the first condition. The cumulative prob- 
ability density P(s) is by definition the quotient of the 
number of normalized spacings Si ^ s by the total num- 
ber of spacings, that is, 



n- 1 



By Eq. (J43J), 

#(s, < a) = #(Ei < £i < ^ £• < £ n ) , (45) 

where 



f± = fi ± v2crWlog 



(46) 



are the roots of the equation s = Wg(£) expressed in 
terms of the maximum normalized spacing 



IT' 



27TCT 



(47) 



Using the first condition to estimate the RHS of Eq. ([45]) 
we easily obtain 

P(s) ~ — [max(f_ -£i,0)+ max(f„ - £+, 0)] . (48) 

In fact, we can replace the latter approximation to P(s) 
by the simpler one 



P(s) 



1 

W 



(£- -£x + £ n - £+) 



(49) 



since the error involved is bounded by 
— \£i+£ n - 2(j\ , 

which is vanishingly small by condition (iv). Substitut- 
ing the explicit expression (|4"6"| for £± into Eq. (|4^)l and 
using (|4^)) and (|4"?| we finally obtain 



P( S ) ~ 1 



The RHS of this remarkable expression depends only on 
the quantity s max , which for the PF chain of BCm type 
is completely determined as a function of N and m by 
Eqs. ©, dlU), dUJ and ([ET ]) - t|B2 |l ■ In particular, for 
large N we have the asymptotic expression 



JhlJl N ^ + o(n-^ 2 ) 

V m + 1 v ; 



Our numerical computations indicate that Eq. (fSTJ)) is in 
excellent agreement with the data for a broad range of 
values of N, m and e = ±1, providing much greater accu- 
racy than the empiric formula (|4"Tj) . For instance, for N = 
20 and m = 2 we have s max = 6^/35/(41tt) ~ 3.12765, 
which differs from the numerically computed maximum 
spacing by 4.5 x 10 -4 . In Fig. [3] we compare the cor- 
responding cumulative spacings distribution with its ap- 
proximation (|50|) using the above value of s max . The 
mean square error is in this case 2.2 x 10 -5 , smaller than 
that of the empiric law (|4"Tj) by more than an order of 
magnitude. 

A natural question in view of these results is to what 
extent the approximation (|50p is applicable to other spin 
chains of HS type. For the PF chain of Am type, one can 
check that the spectrum satisfies conditions (i)-(iv) of 
this section, and in fact we have verified that (f50"|) holds 
with remarkable accuracy in this case [35] . The situation 
is less clear for the originalFIS chain, whose spectrum is 
certainly not equispaced [l6[. Nevertheless, our compu- 
tations show that the formula ([50| still fits the numerical 
data much better than our previous approximation (|4ip , 
a fact clearly deserving further study. As an illustration, 
in Fig. 2] we compare the cumulative spacings distribu- 
tion of the original HS chain with its approximations (|4ip 
and ([50)) for N — 25 and m = 2. It is apparent that the 
new expression (|50p provides a more accurate approxima- 
tion to the numerical data than the empiric formula (I4ip 
(their respective mean square errors are 1.8 x 10~ 5 and 
5.8 x 1(T 4 ). 
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FIG. 3: Cumulative spacings distribution P(s) and its ana- 
lytic approximation ([50} (continuous red line) for TV = 20 and 
m = 2. 




of the Laguerre polynomial lP N 1 as follows: 

(A2) 

The latter sums can be easily computed using the fol- 
lowing identities satisfied by the zeros yi, which can be 
found in Ref. 1331] : 



\ - 2j/j 

E 



Vi - Vj 

( U _ y \ 2 = -Vi + 2(2TV + (3) Vl - (3(P + 4) . 

{Hi Vj) 



(A3) 



(A4) 



Indeed, from the first of these identities we easily obtain 

5> = *(jv+/j-i), E^ = f' ( A5 ) 

i i 

so that, by Eq. ffA4|) . 



E 



2 2/,: 



(yi-yj) 2 6 



= F(-E^ + 2iV ( 2iV + ^ 



W + 4)^i-) = iiV(iV-l) 



Combining the last two equations with (|A1[) and (|A2|) we 
immediately arrive at Eq. (|39[) for the level density /i. 

Turning now to the standard deviation of the energy 
a 2 , from Eqs. (66)-(68) in Ref. [18| we have 



FIG. 4: Cumulative spacings distribution P(s) for the HS 
chain of An type in the case of TV = 25 and m = 2, com- 
pared to its continuous approximations (I50|l (left) and (I41|l 
(right). In both cases, the thin (red) line corresponds to the 
continuous approximation. 



APPENDIX A: COMPUTATION OF /x AND a 2 

In this appendix we shall compute in closed form the 
mean energy \i of the spin chain (TT]) and its standard 
deviation a 2 as functions of the number of particles TV 
and internal degrees of freedom m. 

In the first place, using the formulas for the traces of 
the operators $y , Sij and Sj in Ref. [1 81 ] we easily obtain 



trW e 



iV 



where jj is the parity of m and 



i-^)5>> (ai) 

771/ ' 



The sums appearing in Eq. (|A1|) can be expressed in 
terms of sums involving the zeros yi = £^/2, 1 i TV, 



2_trH? 2 



?7i 



N 



i-^)(2E(4+^)+E^) 



Ad-rf(iE^-E^A 



(A6) 



All of the sums appearing in the latter expression can be 
readily evaluated. Indeed, we have 



4 ^ (yi-yj) 2 i6(/? + i) 



(A7) 
(A8) 



where we have used Eqs. (15) and (17) from Ref. [3c 
On the other hand, 

<*y <*y 12/1 %j 



y? + y 3 2 3 x - l 



2 E 



4 E 



//r 



3TV(TV- 1) 
(y.-j/,) 4 " 8(/3+l) ' 



(A9) 
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while from [34|j Thm. 5.1] it follows that 

+ 2N(8N 2 + 2(4N - 1)(3 + 3/3 2 ) 
- 4(2N + f3)((3 2 - 2/3 - 18) V — 

+ f3((3 3 - 4/3 2 - 104/3 - 144) ^ . (A10) 

All the sums in the right-hand side of the latter expres- 
sion have already been evaluated in Eqs. (|A5[) and (IA7|) . 
except the first one. In order to compute this sum, we 
multiply Eq. (|A3|) by yi and sum over i, obtaining 



E 



2yf 



E 



2 2 



= E^ i + %) 



= 2(JV -!)$> = 



and hence, by Eq. (|A5|) . 

^y?=iV(iV + /3-l)(2iV + /3-2). (All) 



Substituting the value of the latter sum in Eq. (|A10[) and 
using (IA5D and (IA7D we obtain 



AT (AT - 1)((27V + 5)/3 + 2A + 14) 



\ - yt 



144(/3+l) 



(A12) 



Equation (|4T))) now follows by inserting (|A7|) . (|A8 
and fX9 ]) -(rXT2 | into Eq. (|A"6|) . 



APPENDIX B: COMPUTATION OF THE MINIMUM 
ENERGY 

In this appendix we shall obtain an explicit expression 
for the minimum energy S^in °f the spin chain JTJ). Our 
starting point is Eq. ([7]) , which implies that £^ lin is given 
by 

^min = 1™ — (-E m ; n — E min J 

a^oo a 

in terms of the minimum energies E^ in and E^ in of the 
scalar and spin dynamical models ([3]) and (J5]), respec- 
tively. By the discussion in Section [II] (cf. Eq. (fT§)) ). £^ lin 
is the minimum value of |n|, where n is any multiindcx 
compatible with conditions (i)-(iii) in the latter section. 
From these conditions it follows that the multiindex n 
minimizing |n| is given by 



: k : k 1 . 



,o,...,o) 



where ^ I < Ik and U is given by 
- , m even 



k 



odd and i even 



"^Y^ j m an( i i odd . 



In view of the above expression, it is convenient to treat 
separately the cases of even and odd to. For even to, we 
have li = ^ for all i, so that 



k=[2N/m\, l = Nmod—, 



and 



fe-i 



^i + lk = — k(k- 1) + lk 



N 2 N l(m-2l) 



to 2 



2m 



(m even) . (Bl) 



Suppose now that to is odd. If k = 2j is an even number, 
then l = l 2 = --- = l 2 j-2 = h = h = ■ • • = hj-x = 



- and thus N = jm + I, so that 
j = [N/m\ , I — (N mod m) < hj 



m + e 



The minimum energy in this case is thus given by 

3-1 3-1 
to + e sr~* to — ev-i, 

£min = —j- E 2 * + E( 2i + !) + 2 -? 7 

i=0 t=0 

to — e 

= - 1) + J — 5 h J 

_ N (m + e) Z(m + e-2Z) 

to 2m 2m 

On the other hand, if k = 2j + 1 is odd then Zq = Z2 



Z 2 j — 1 = 2 V £ and thus 



; rn+t ; ; 

••' — hj — 2 ' '1 — f 3 — • • 

Ar = jm + I + 2S±£, with < Z < l 2j +i = ZfZ. Callin 
Z' = I + ^ we have 



j = [N/m\ , I' = (A" mod to) ^ 



TO + £ 



and 



. 3 3-1 

m + e v-^ m — ev-^, 

^ni„ = -5- E 21 + E( 2i + !) 



i=0 



+ (2j + !)(/'■ 



i=0 

to + e 



TO — 6 , , TO + e 

= mj(7 - 1) + j + 2jl' + 1' — 

_ N 2 N(m + e) (I- m)(m + e - 21) 

m 2m 2m 

Hence we can express the minimum energy for odd m in 
a unified way as 



£ e . = 



N 2 N(m + e) 1 



+ — (l-m8(2l 



0) 



to 2m 2to 

x(m + e-2Z) (modd), (B2) 
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where 

I = N mod m 
and 6{x) = 1 for x ^ and otherwise. 
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